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CHAPTER I 

introduction 

To meet the performance demands of future military 
and commercial aircraft, propulsion system requirements 
will make necessary the design and development of sub- 
stantially more complex engine configurations. These new 
propulsion systems will dictate the need for advanced multi- 
variable control system design methodologies capable of 
interrelating the design goals with ths" physical limitations 
of the turbomachinery. This means that although steady 
state operating line conditions will retain their importance 
in the design specifications further emphasis will be placed 
on the short term transient performance of the engine. To 
satisfy these performance requirements, more control vari- 
ables will be introduced and more output variables will be 
monitored to evaluate engine performance. 

in the early years of jet aircraft propulsion, single- 
spool turbojets were used almost exclusively. The first 
control systems for these engines wre hydromechanical and 
used the principle of the flyball governor extensively 


for fuel-rotor speed control [1] . As tubojet engine 
performance demands increased, a second compressor- 
turbine section was added to drive a fan to remove energy 
from the primary jet airstream* thus reducing engine air- 
flow velocity and increasing cycle efficiency. The addition 
of this second compressor-turbine section enabled greater 
flexibility of compressor performance at high discharge 
to intake pressure ratios, but required the control of two 
mechanically independent turboshafts . Turbofans now 
power most modern subsonic aircraft and, with afterburners, 
many high performance supersonic military aircrafts as well. 

The turbofan engines in use today are basically fixed- 
geometry designs. The control system modulates one or two 
variables and is configured as a set of single loop controllers 
which are scheduled in an open loop manner according to 
aerodynamic and environmental conditions. The majority of 
the control systems have been designed using single-input 
single-output (SISO) classical frequency response techniques . 
The resulting control configurations have typically been 
implemented using hydromechanical or electronic analog 
components. 

As aircraft turbine engines become more complex, 
additional requirements for the control system will result. 
Operational requirements will become more stringent thus 


/ 




necessitating increased control system capability in terms 
of computational accuracy and static and dynamic response. 

As performance requirements increase, more variables must 
be controlled and more engine parameters must be monitored. 

SISO control system design methods will no longer be 
appropriate as cross-coupled system interactions must be 
taken into consideration. Open-loop, trajectory-dependent 
control designs will yield to closed loop design procedures 
as the now-standard, analog/hydromechanical devices yield 
to full authority digital electronic control packages. 

Some research and development activity for turbofan 
engine control systems using modern control theory has been 
reported by Michael and Farrar [2,3] , Stone and co-workers 
[4], Weinberg [5], and Merrill and Leininger [6]. The most 
comprehensive application of modern control theory to air- 
craft propulsion system design was performed by Systems Control 
Inc . on the Pratt and Whitney FlOO-PW-lOO turbofan engine [7- 
11] . This LQR control design was successfully tested in an 
altitude chamber at NASA-LeRC by Lehtinen and co-workers 
[12,13] . The results of this test demonstrate the effectiveness 
of modern control theory for highly non-linear sophi- 
sticated turbofan applications in addition to providing 
an evaluation of digital electronic devices for turbofan 
engine control systems. 


In 1977, the FIDO sea level static condition was used 
as the theme problen\ to exploro alternatives for linear 
multivariable control in a symposium sponsored by the National 
Engineering Consortium. The results reported in [14] demon- 
strate the successful application of multivariable frequency 
domain control system design methods to the PlOO design 
problem. 

During the period 1974-1979 , the General Electric 
corporation, under contract to NASA, designed, developed, 
and tested an experimental commercial aircraft engine [15,16] . 
This new Quiet Clean Shorthaul Experimental Engine (QCSEE) 
incorporates performance and structural characteristics 
unlike those in any turbofan engine in production today. 

The QCSEE program introduced high bypass ratios for turbo- 
machinery noise suppression, reversible pitch fan blades for 
rapid thrust response, total digital electronic engine controls, 
and extensive use of composites for drag reduction and weight 
considerations* In addition to fuel flow and fan blade pitch 
angle, variable nozzle area control is also provided . To 
incorporate all of these characteristics into a single 
propulsion system, represents a significant breakthrough 
in turbofan engine technology. 

During the period of development and testing of the 
QCSEE engine, NASA developed a highly non-linear, accurate. 







real-time digital simulation of the engine at sea-level 
static conditions. This non-linear model was used in ex- 
tensive tests at the NASA-Ames in-flight simulator facility 
for test pilot evaluations of integrated engine airframe 
combinations (171. 

Using the non-linear digital simulation as a repre- 
sentative model of the dynamic operation of the QCSEE 
turbofan engine? a feedback control system is designed in 
this report by multivariable freguency domain design tech- 
niques. Transfer function matrices are generated at each 
of five power lever settings covering the range of operation 
from approach power to full throttle, i.e. , 62.5% to 100% 
full power. These transfer functions are then used by an 
interactive control system design synthesis program to 
provide a closed loop feedback control using the multivariable 
Nyquist array (MNA) method [18-20] and the recent extensions 
to multivariable Bode diagrams (MBD) and Nichols Charts 
(MNC) [21-22] . 

The QCSEE design is initially performed holding nozzle 
area full open with fixed fan blade pitch angle at the power 
setting of 62.5% of full power. This SISO design is evalu- 
ated at other power settings and tested in the non-linear 
simulation to evaluate engine performance during a power 
slam from 62.5% to full power (100%) . 


The second phase of the MNA design uses fuel flow and 
fan blade pitch angle in f ealizing a two input two output 
feedback design. These results are extended to the three 
input three output system with nozzle area as the third 
input. This variable is used to provide additional control 
over inlet~Mach number for noise suppression. In each case^ 
the MNA design is compared with the Standard QCSEE control . 

The next chapter details the basic concepts associated 
with the Multivariable Nyguist Array Method and introduces 
several recent extensions to the theory. Chapter 3 then 
uses this design philosophy to determine the closed loop 
feedback control system for the QCSEE turbofan engine. The 
successful application of the MNA method clearly demonstrates 
the effectiveness and utility of the multivariable classi- 
cal frequency design procedure. 
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CHAPTER II 

MULTIVARIABLE NYQUIST ARRAY METHOD 

The Multivariable Nyquist Ai'ray (MNA) method is com- 
posed of two Nyquist array design procedures : one pro- 

cedure for the direct polar plane (DNA) and one for the 
inverse polar plane (INA) . Both methods have identical 
design objectives and are founded upon a common mathematical 
structure. The principal point of departure is the use and 
interpretation of the multivariable Nyquist stability 
criterion in achieving the final system design. 

The INA method was first introduced in 1968 [23] based 
upon the premise that compensator selection to achieve the 
desired diagonal dominance condition was easier to obtain in 
the inverse polar plane. In 1974, a more comprehensive 
description of the INA method became available [24] which 
also provided the forum for the introduction of the DNA method. 
Until recently these design methods were considered to be 
mutually exclusive in the sense that once the open loop 
design conditions were met in a specific polar plane the 
closed loop design analysis was restricted to that polar plane. 
Thus an INA design could not at that time be projected into 
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tho more favorod direct polar plane, fhin led Fisher et al 
[25] in a comprehensive application study to indicate a 
preferencQ for the DNA method over the INA '.Method whereas 
the INA method is used most frequently in Europe. 

In this chapter the MNA design objective is briefly 
reviewed. Now theoretical concepts are presented which 
eliminate the polar plane dependency and establish a found- 
ation for finite frequency range designs. A new graphical 
design procedure is then introduced which utilizes the 
fundamental ideas of Bode and Nichols in multi-loop systems, 

A. MNA Background and Notation 

For the system in Figure la, G(s) represents an mxm 
transfer matrix relating the output vector y(s) to the system 
input vector u (s) . The vector y^.^^ (s) denotes the reference 
output vector. The design objective is to determine the input 
vector u(s) such that the system output vector follows the 
reference vector with dynamics acceptable to the system 
designer. If, for example, yj-^fCs) represents a desired set 
point Operation then y (t) should achieve the desired level 
within prespecified rise time, settling time, and percent 
overshoot requirements. 



To the ssystem of Figure la, three square compehsators 
are used to close the feedback loop as in Pigure lb. For 
reasons which will be apparent later in the development 
li and P(s) are assumed to be diagonal and K is assumed to be 
an arbitrary constant matrix, All matrices are assumed to 
be of rank m. Shifting the constant diagonal matrix h to 
the loft aide of the summing junction Figure Ic is obtained. 
Re-arranging Figure Ic and letting 


^rof “ ^ ^ref 
Z(t) » h y(t) 



the standard MNA block diagram of Pigure Id is obtained. 
Let P(s) equal the identity matrix and define 


Q(s) * L G(s)K (2) 

Following Rosenbrock [24) , K and L must be selected such 
that Q(s) and/or 

Q(s) = K G(s)L (3) 

is diagonal dominant by rows or columns for all s on the 
standard Nyquist D contour. 



















The concept of diagonal dominant matrices is well 
documented [18,23r24] and is defined as follows for the 
DNA method: 


Definition i 


Let Q(s) be a complex square matrix 
of dimension m. Then Q(s) is 
diagonal dominant if either or both 






for 

all 


m 

= Max 

I 

seD 

j=l 




m 

= Max 

I 

seD 

j-1 




(s) I < 1.0 


Iqji(s) |/|qii(s) I <1.0 


(4) 


(5) 


where (4) utilizes the rows of 
Q(s) (row dominance) and (5) 
uses the columns of Q(s) (column 
dominance) . 


Clearly if Q(s) in (2) (or Q(s) in (3)) is purely 
diagonal, then m feedback controllers (and m loop gains) 
can be designed independently for the uncoupled system. 


To achieve this di?gree of independence, K and L would of 


necessity be complicated functions of s and could result 
in severe implementation problems. Rosenbrock guestioned 
the need for a purely diagonal open loop transfer matrix Q(s) 
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and developed a procedure wherein the off-diagonal entries 
In Q(s) were allowed to be non-zero but of sufficiently 
small moduli to permit design synthesis by single loop 
methods. By relaxing the diagonality condition of Q(s)i 
to that of a diagonal dominant 0(s)f the compensator matrices 
K and L in most cases remain constant thus significantly 
improving the implementation requirements. 

Using the above definition for diagonal dominance ^ 
Rosembrock was able to extend the classical SISO Nyquist 
Criterion to the multivariable case. This development led 
to a graphical interpretation of the multi-input multi-output 
(MIMO) Nyquist design objective where the traditional 
Nyquist curve is replaced by a band of circles. The radii 

i 

of the circles (called the Gershgorin circles) consist 
of the sum of the moduli of the off-diagonal terms in Q(s) 

by rows or by columns. Pot the DNA method the Gershgorin 

radii for ( 4 ) and (5) become 

■ m 

= I ki-: (s) I (6) 

j=l ■' 
jj^i 

_ m 

= I lq..(s)| (7) 

j=l J 















13- 


Henoe, the dominance levels can be written as 


0^ = Max d|(s)/lq. . (s) | 

(8) 

scP 


0? = Max d1(s)/|q^^(s) 1 

(9) 


To determine the elements of K and L, an optimization 
procedure was developed by Leininger [18,19] based upon 
the dominance levels in (8) and (9) » Pundamentally, the 


optiralzation objective is to secure the lowest level of 


dominance in the respective rows (or columns) using a 
constant compensator structure. This method of achieving 
dominance can be used in conjunction with any of the twelve 
MNA forms and thus significantly reduces the time required 
to obtain dominance from several weeks (or months) using 
previous trial and error procedures [24] to several CPU 
seconds. 

Once the compensator pairs have been selected in 
accordance with the above concepts, the respective control 
loops can then be designed independently using classical' 

SISO design procedures. The resulting closed loop design 

configuration can be further reviewed on ah individual 

loop basis through the utilization of Ostrowski' s theorem [26] 



, This theorem provides for the reduction of the Gershgorin 

radii in loop i based upon the feedback configurations 
in loops j^i* Thus a more accurate assessment of the respective 
I phase margins and gain margins cm be made. 

It is useful at this point in the development to dl- 
l' gress from the formalities of the fINA design procedure to 

I closely examine the closed loop structure that has been 

I 

I created above. This analysis will place the MNA design 

1 procedure in proper perspective and provide the foundation 

for subsequent designs using multivariable Bode and Nichols 
I methods. 

I 

B. MNA - Perspectives 

In most applications of the MNA design method Q(s) 

(or Q(s)) cannot be made diagonal dominant when L is the 
identity matrix and K is constrained to be a constant 
matrix. The alternatives are therefore to eith®^ insert 
a constant postcompensator matrix or allow K to be a dy- 
namic function of s. To design a K(s) for open loop 
dominance over a large frequency range would be an ex- 
tremely cumbersome and tedious task. It is also questionable 
whether the end result (i.e, dominance of Q(s) ) could 
justify the means by which it was achieved . That is, if 


(■ 
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dynamics are to be Introduced into the feedforward path, 
it should be to the benefit of the closed loop design 
rather than achieving a particular open loop conf iguratipn. 

The problem of inserting loop dynamics to meet design 
specifications will be addressed in the next subsection^ 
Consider the utilization of a postcompensator matrix 
L as indicated in Figure 1-C# Note that the shift of matrix 
L to the left side of the summing junction completely 
changes the characteristics of the design problem unless 
L is constrained to be diagonal. For a non-diagonal L 
the design problem becomes one in which a linear combination 
of the system outputs, L y(t), are expected to track a linear 
combination of the reference values, h yj.gj. The error 
vector for this case is thus a linear combination of the 
output errors , 

Oy(t) = L(y^g£ - y(t)) (10) 

and may become zero for conditions other than 

y(t)=>y^gf (11) 

However , when L is diagonal the error vector in (10) becomes 
zero if and only if (11) is satisfied. We therefore, im- 
pose this constraint on L and utilize the concept of dominance 
sharing [19,20] to achieve the desired dominance condition 
for Q(s) (or Q(s) ) . 
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In Figure Id^ assume F(s) is the identity matrix^ 

L is diagonal/ and K and L are constant matrices which 
have been determined such that Q(s) in (2) is column 
diagonal dominant * A common misconception associated with 
MNA methods concerns an alleged modification of the rela- 
tionship between y(t) and u(t). Clearlv / this relationship 
is always dictated by G(s) and remains intact regardless 
of the selection of K and L. What is dependent upon K and 
h, however/ is the relationship between the corresponding 
entries in z(t) and e^ (t) , 

2 (s) = Q(s) (s) (12) 

When Q(s) is diagonal dominant, 2 j^(s) is primarily dependent 
upon the i th error signal 

= 'ref - <“) 

with cross"Coupled error signal (ij^j) interactions determined 
by the degree of dominance achieved and respective loop gains 
(inherent within the elements of Q(s)). 

The i th element in (12) can be written as 

m 

z.Cs) =q..(s) e (s) + I q..(s) e (s) 

X XX j=i. Zj 

j?^i 


( 14 ) 
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where i » 1, 2 , . . . , nii Noi:e that: although we assumed 
Q(s) column dominant, it is the row elements of Q(s) which 
contribute to the interaction terms in (14). To examine 
this anomaly further, let m “ 2 and postmultiply (2) by 
P as in Figure Id 

Q(s)P « L G(s)K F (15) 


where F is a diagonal constant gain matrix. Since Q(s) 
is assumed column dominant, the dominance levels of Q(s) 
are unaffected by the moduli of the elements of F. Thus 
(12) becomes 


Z3_(s) - gj^l (s) f^ e^^(s) + ^2 ®Z2 

Z2(s) = q2i(s) e^^(s) + q22{s) f2 (17) 


With lq..(s)| > jqj.(s)l for all seD it is clear that the 

XX J X 

error term e„ (s) has proportionally more influence on z.(s) 

'‘i ^ 

than on Zj(s), ij^j . But note what happens as the gain f 2 

is increased with f, fixed. The interaction of e (s) on 

* ' “2 
z. (s) increases while the effect of e_ (s) on z_(s) decreases. 

X ^ 

This interaction influence of the gain matrix F is independent 
of the respective dominance levels. 

Thus low levels of dominance associated with each 
diagonal entry in Q(s) are not sufficient to insure low inter- 


action levels unless the impact of the respective loop 
gains has been properly addressed. This analysis explains 
why the DNA (or INA) bands may be very narrow (lov? dominance 
levels) yet the closed loop time domain responses indicate 
high levels of cross-coupled interactions. 

To pursue this analysis further, consider the closed 
loop system of Figure Id and assume that K and L have been 
determined for column dominance of Q(s). Shift the diagonal 
L matrix around the loop as in Figure 2 recognizing that this 
is the most probable implementation configuration since 
our interest is in the y(t) vector and not necessarily z(t) 
where 


z(t) = L y(t) 


(18) 


From Figure Id, the closed loop relationship is 


-1 


Z(s) = (I + LGKF) LGKF Zj.^^(s) 


(19) 


or 


Z(s) = P(s) Zref<s) 


( 20 ) 


with 


P(s) = (I + LGKF)"^ LGKF 

= (I + Q(s)F)"^ Q(S)F 


( 21 ) 





( 22 ) 
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Using Figure 2 the closed loop system is described by 
y(s) (I + G{8) KFL)“^ G(s) KFIi y,.ef<s) 
or 

y(s) » H(s) 


H(s) = (I + l“^ Q(s) FL)“^ Q(s) FL 
Comparing (21) with (24) if follows that 


or 


H(s) - l"^ P(s) h 


H(s) = 


Plj^(s) 


Poi (S) 


Z- ^ 21 


^12^®^ • * 


^22 ^®^ 




I~ ^ml^®^ 
in 




^<®> 


(23) 

(24) 

(25) 


(26) 


Hence, the diagonal elements of the closed loop transfer 
function matrix are unaffected by the compensator shift. 
Further, if each control loop is designed for a small 
steady state error condition (i.e. , high loop gains), 
then the off-diagonal terms in ( 26 ) tend to zero and y(t) 
tracks If the gain in loop i decreases , then the 

degree of interaction in loop i from the remaining high 
gain loops will increase. This interaction effect is not 
easily detected from the Nyquist array diagrams as it is 
independent of the dominance levels when using a DNA * 
column or an INA row dominant design objective. 

The concerns associated with closed loop interaction 
levels and system dominance can be resolved simultaneously 

when using a DNA row and/or an INA column dominance objective. 

2 

For these MNA design configuration the m elements of K, 
m elements of L, and the m dynamic compensators of F(s) 
must be simultaneously selected. In this way, system domi- 
nance, loop gains, and response shaping for all loops are 
addressed simultaneously since they are highly interconnected. 
As a result of the complexity of these design formats, 
they have never been seriously considered as viable design 
procedures. 

In the following section recent extensions to the MNA 
design technique are briefly described. 
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C. New Design Procedures 

The MNA method as originally conceived by Rosenbrock 
was composed of two design procedures. Each design pro- 

f 

cedure was relegated to the specific complex plane in 
Which the diagonal dominance conditions was obtained. 
Stability theorems based upon the extended Nyquist cri- 
terion were then developed for each design procedure. » 

In a typical application of the MNA design methods, 

Q(s) and Q(s) are not, in general, simultaneously diagonal 
dominant in the Gershgorin sense for all s on the Nyquist 
D contour. This is due to the fact that 

(s) / qjj (s) (27) 

for i, j , and s on D. Hence it was thought that the MNA 
forms were not compatible, i.e. , DNA compensation methods 
could not be used when INA dominance was obtained and visa 
versa. This limitation of the MNA design philosophies 
is removed by introducing a conformally mapped image process. 

The details of the imaging process are contained in 
[20] and will be briefly reviewed . 

Assume Q(s) is Gershgorin dominant in the inverse 
polar plane Z where the origin is exterior to the respective 
Gershgorin bands. Let 


W = 1/Z 


(28) 
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represenib the conformal mapping from the inverse polar 
plane to the direct polar plane W. Thus under (28) there 
exists a one-to-one correspondence between the points 
z(s) * x(s) + j y(s) in the Z plane and the points w(a) 
u(s) + j y(s) in the W plane. Since (28) is a conformal 
mapping circles in the Z plane map to circles in the W 
plane. 

Let d(s) represent the radii of the cirlces in the Z 
plane and r(s) represent the radii of the image circles 
in the W plane. Then the image center point and image 
radii are given by 


u(s) 


x(s) 

x^(s) + y^(s) - d^(s) 


v(s) a 


x'^Cs) f y'^(s) - d'^(s) 


d(s) 


x^(s) + y^(s) - d^(s) 


Since the origin is not an interior point of the Z plane 
circles, the interior of these circles map to the interior 
of the image circles. Furthermore, the relationship be- 
tween respective circles in the W plane will be maintained 


Bince (28) is conformaX, Hence, any gain and/or phase 
margins measured with respect to the Gershgorin circles 
(or envelope) in the 2 plane will be retained in the image 
plane. 

'This mapping procedure works egually as well when 
Q(s) is diagonal dominant and the Gershgorin bands are 
projected into the inverse plane. Note that the imaging 
process is independent of Gershgorin dominance levels in 
the image mode as the imaging process retains the dominance 
levels from the projected domain. ' 

Stability theorems for the image bands as well as 
further insights into the imaging process are contained in 
Appendix A. The significance of this theoretical development 
will become more apparent when the Bode and Nichols methods 
are discussed for MIMO systems later in this section . 

A second limitation of the MNA design methods as pro- 
posed by Rosenbrock was the requirement of diagonal dominance 
for Q(s) or Q(s) by rows or by columns over the entire 
Nyquist D contour . This condition was imposed to ensure that 
the stability theorems of Rosenbrock would remain valid. For 
if the system lost diagonal dominance, then the theorems 
neither confirm nor deny stability. 
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Since all jreal systems are band limited (i.e., strictly 
proper in Rosenbtock's terms) , the search for dominance 
only needs to be checked as w proceeds up the imaginary 
axis on the D contour. Although this simplifies the search 
to the customary range of 0 £ w £ ~ it imposes an excessive 
burden on the system designer if he knows that frequencies 
beyond some finite value, are not likely to contribute 
to the system response. 

In response to this concern, a new theoretical develop- 
ment was introduced wherein the contribution to the MNA or 
image bands for all frequencies beyond w^ could be properly 
accounted for. The specific mechanism for finite frequency 
domain design eval'^at ion and related stability theorems are 
detailed in [20] . 

Fundamentally, the contribution to the dynamic response 
for frequencies beyond some w^ are contained within a circle 
centered about the origin in the complex plane. The radius 
of this circle is determined by the moduli of the elements 
in Q(s) or Q(s) at s - w^ increases the corresponding 

circle radii decrease for the DNA method and increase for the 
INA method. The feedback gains are then selected in accordance 
with the appropriate array bands and origin centered circles. 


The origin -centered circles associated with a finite 
frequency range evaluation can be conformally mapped to 
the image plane. In this case the origin is contained within 
the circle, hence under the mapping principal the interior 
of the circle in one domain maps to the exterior of the 
circle in the image domain. ^ 

The two new theoretical developments discussed above 
provide the basis for the equivalence between the multi- 
variable i^yquist array design procedures. As a result, 

i 

closed loop compensator design analysis and synthesis is 
no longer restricted to the polar plane in which the desired 
diagonal dominance condition was obtained. Thus, DNA 
designs can be used for INA dominance conditions and the 
desirable INA minor loop synthesis procedures can be used 
when DNA dominance has been achieved. In addition the 
theory is now complete for finite frequency dominance, 

Nith this equivalence principal we are now in a position 
to extend the design synthesis procedure to the ma!gnitude/ 
frequency and phase/frequency planes. The resulting graphi- 
cal display will allow for design analysis and compensator 
synthesis using Standard Bode techniques. Using the 
raagnitude/phase angle piane, a graphical display for each 
control loop can be reviewed by the familiar Nichols chart 
methods . 


INA DIAGRAM I BODE-DB 



Figure 3: Typical INA- IMAGE - BODE-MNA PLOTS 


The fundamentaX concept of multivariable Bode Plots 
and Nichols Charts to be presented below are described in 
more detail in [21] . Reference [22] describes a computer 
aided design suite for closed loop system design and includes 
all of the graphical procedures referred to in this chapter. 

The multivariable Bode diagrams are constructed 
directly from the data used to calculate a direct polar 
plane Nyquist diagram. If the OLTF matrix is DNA dominant, 
the envelope curves [27] of the corresponding Gershgorin 
radii are each converted to polar form and independently 
drafted as a function of frequency. When system dominance 
is obtained in the inverse polar plane, it is necessary 
to first form the equivalent set of image bands in the 
direct polar plane. The envelope of the image bands are 
then used to construct the multivariable Bode diagrams. 

Figure 3a represents an INA diagram which is dominant 
over the finite frequency range 0 £ w £ 30. Figure 3b 
contains the corresponding image band for Figure 3a with 
Figures 3c and 3d displaying the Bode information. 

Since the mapping between the respective complex polar 
planes is confojnnal, it is not necessary to actually compute 
and/or display the image circles. The envelope curves for 
the image domain can be computed directly from the envelope 
curves in the domain in which dominance was achieved . 


Theorem 


Let and represent the envelope curves 

for an MNA diagram. The envelope curves, C.(w) and C 2 (w), 
to the corresponding image circles in the image plane are 
given by, 

C^(w) = tCj^Cw)]"^ and C 2 CW) - [C 2 (w) ] (32) 

proof ; 

The relationships in (32) follow directly from the 
conformality of th® ™apping W = 1/Z between the respective 
polar planes. 

Q.E.D. 

Thus for an INA system dominant condition, the Bode 
diagrams can be calculated based upon the inverse of the 
INA Gershgorin bands. The multivariable Nichols chart 
is then obtained directly from the data used to construct 
the Bode diagrams and is presented in Figure 4 for the 
conditions detailed in Figure 3. 

With the graphical analysis procedures described above, 
system output response shaping can be considered by any of 
the standard SI SO methods. Since INA and DNA design methods 
are described elsewhere, the remainder of this section will 
be restricted to the Bode design form. 





Consider a DNA column dominant condition for Q(s) 
in (2) , and assume the system is type o. Since the 
loop gain determined at w = 0 may be low, a reasohably 
large steady state error could be expected. If the 
design specification will accept an error no greater 
than 10% , we may compute the required gain increase 
using 

®ss ** 1 + f^ P^(0) (33) 

where P^(0) is the value of the i th Gershgorin envelope 
(or image DNA envelope) closest to the origin at w = 0. 
Using this value in (33) ensures that the steady state 
error will always be less than the required minimum level. 
Thus 


fj^ = [— ^--^]Pi(0^ i = l, 2, ..., m (34) 

\ ss 

Similar relationships can be obtained for type 1 or 
higher systems. 

Inserting (34) into the system will only shift the 
magnitude curve of the Bode diagram. As the curve is 
shifted, the closed loop bandwidth of the system (-3 db 
point on magnitude curve) will be simultaneously affected. 
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It is also possible the system dominance may be lost 
as a result of the gain change. This is apparent for 
those cases in which the 180“ phase angle information 
corresponding to the new sero db crossover point is 
contained within the phase envelope curves. 

This condition, howevetj does not present any 
serious complications since we can easily synthesize an 
appropriate compensator to shift the -180“ point outside 
the phase angle bands at the crossover frequency* This 
can be accomplished using standard lead, lag, or lead-lag 
methods. The details of the compensation procedure 
are described in t21, 22], Similar procedures apply for 
type n systems. 

It is now apparent that the multivariable Bode diagrams 
and Nichols chart can play an important role in the effective 
design of multivariable closed loop systems. In addition, 
it provides an effective measure of the frequency range 
over which dominance was obtained. If the desired closed 
loop bandwidth is beyond the range in which dominance was 
obtained, new efforts are required to compensate for this 
deficiency. At the other end of the spectrum, if the zero 
db crossover frequency is well within the dynamic range, 
this could reflect too conservative of a frequency range 


(i.e., frequency range too large) for dominance acquisi- 
tion. New design efforts over a reduced range may provide 
less interaction levels in the bandwidth regions of interest 
and could improve the overall design significantly. 

D, Summary 

In this chapter the MNA design procedure was briefly 
reviewed and new design tactics were introduced. The Image 
bands play an important role in establishing the equivalence 
of the INA and DNA design methods and provide the theoretical 
foundation for viewing the INA Gershgorin envelope curves 
in the Bode and Nichols domain. The finite frequency de- 
velopment provides the theory needed to reduce the com- 
putational burden for the diagonal dominance search. The 
multivariable Bode diagrams provide a new view of the system 
design where SISO compensator design procedures can be 
used effectively to secure dominance and simultaneously 
achieve the desired closed loop performance. The correspond- 
ing multivariable Nichols chart displays the Bode information 
in a familiar form which again is suitable for design 
analysis and control system synthesis. 

The next chapter utilizes the design procedures and 
theoretical developments of this chapter to design a closed 
loop control for the General Electric QCSEE turbofan engine. 


CHAPTER III 


MNA APPLIED TO QCSBE ENGINE 

The turbofan engine selected for the MNA control 
design study was the quiet clean shorthaul experimental 
engine (QCSEE) [16] . The QCSEE engine program was a NASA 
sponsored effort to design, build, and test two experimental 
engines which would serve as test beds for advancsd 
propulsion technology evaluation. One engine was designed 
for Under “the-wing (UTW) mounting while the second design 
was developed for an over-^the-v^ing (OTW) mounting. This 
study was limited to the UTW design. 

Turbofan engines of the QCSEE design are under con- 
sideration to power shorthaul aircraft of the STOL con- 
figuration. These aircraft would operate from small air- 
ports located near major urban communities. In this way 
large transports could land at metropolitan airports mhny 
miles from the urban center with passengers and cargo 
transported to the civic centers via low noise and high 
reliability STOL aircraft. 

To provide these services, primary propulsion re- 
quirements for the aircraft are those of low noise levels 
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and exhaust emissions at all operating levels, rapid 
thrust response capability including reverse thrust and 
low weight. The design challenges were met by the 
General Electric Company. 

In the next section the specific design for the OCSEE 
UTW engine is detailed with principal emphasis placed on 
the control system design requirements. The following 
section describes the non“linear digital simulation program 
developed by NASA from which the linear operating point 
models for the MNA design are obtained. The last section 
discusses the results of the control system design study 
performed using the multivariable frequency domain design 
procedure. 

A. QCSEE Turbofan Engine 

A cross-section of the under-the-wing QCSEE turbofan 
engine design is detailed in Figure 5. The design in- 
corporates performance and structural characteristics 
unlike those in any engine in production today and includes 

1. An extremely high by-pass ratio (12:1) 
and a high throat Mach number inlet for 
noise suppression. 

2. Reversible pitch fan blades for rapid 
thrust response. 


QCSEE UTW ENGINE 



3. Geared turbine/fan combination for low 
fan tip speeds with a high thrust rating. 

4. Digital electronic engine controls. 

5. Extensive use of composites for drag 
reduction and weight con'siderations. 

The engine also incorporates a variable area fan duct 
exhaust and variable compressor stator blades. The fan 
is driven by the low pressure turbine through reduction 
gears to achieve an efficient low pressure ratio# low 
tip speed operation. The variable pitch mechanism is 
capable of moving the blades from a forward thrust position 
to a full reverse thrust position in less than one second. 
In the reverse thrust position, the hydraulically actuated 
fan duct exhaust nozzle is opened to a flare position 
to serve as an efficient inlet. 

The power requirements of a shorthaul aircraft system 
dictated the design specifications which had to be met by 
the QCSEE control. The primary specifications were: 

1. to allow the pilot to command percent 
rated thrust, 

2. to assure engine operation at all times 
within safe engine limits, 

3. to control inlet Mach number to meet 
noise suppression levels, 

4. to achieve rapid thrust response (in 
both forward and reverse modes). 



These requirements all relate mainly to safety and noise 
emission considerations. The first requirement resulted 
from the desire to decrease pilot work load in a STOlj 
environment. This can be achieved by closing the feedback 
loops using a parameter closely related to thrust (such 
as engine pressure ratio) instead of using a parameter 
indirectly related to thrust (such as spool speed) . This 
reduces the time the pilot must spend adjusting the 
throttle to achieve the desired thrust level. 

The second requirement meant that logic must be 
incorporated to prevent compressor and fan overspends 
and turbine overtemperabure . To meet noise suppression 
specifications, a high inlet throat Mach number is required 
to significantly reduce noise radiation from the inlet 
duct. The requirement for rapid thrust response arises 
from the need to insure safe STOL operation. In particular, 
to have a safe go-around capability, thrust changes from 
62% (approach power) to 95% of full power must be achieved 
within one second. This is two to four times as fast as 
is possible with current production engines. 

The above control system requirements were met by the 
engine manufacturer by using an engine mounted digital 
control system shown in Figure 5. All primary control 
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functions are progranuned on the digital control computer. 
Digital command signals are sent directly to the three 
input actuators: 

1. fuel flow actuator, 

2. fan pitch actuator, 

3. nozzle area actuator. 


A hydromechanical control provides fuel flow limiting and 
is also used to provide the command to the compressor 
variable stator actuators. 

An important consideration in the manufacturer’s 
design process was the control mode selection. This is 
the process of pairing the four control variables with 
four output variables. The mode selection was made such 
that at key, selected steady-state operating conditions, 
primary (unmeasurable) performance variables such as thrust, 
specific fuel cohsumption, and stall margins are relatively 
insensitive to small sensor inaccuracies. The result of 
the mode selection process for the QCSEE UTW engine control 
was to pair fan speed with fan pitch angle, engine pressure 
ratio with fuel flow and inlet Mach number with 

fan duct nozzle area. The core stators were controlled 
essentially independently, being scheduled (according to 
standard practice) as a function of corrected core compressor 
speed . 


1 


The steady state operation of the QCSEE engine is 
determined by a set of demand schedules. These schedules 
command optimum values of engine pressure ratio (EPR ) , 
inlet Mach number and fan speed as functions of the pilot's 
power request and ambient conditions. In particular, the 
EPR schedule is a monotonic increasing function of the 
power request. The inlet Mach number is scheduled to be 
a constant 0.79; however, at low power (less than 701 of 
maximum) , the nozzle is full open and the Mach number then 
varies with power request. The fan speed is scheduled 
as a constant value for power requests above 55%. Below 
that value, fan pitch is at its minimum allowable setting 
and hence fan speed varies with power request. Keeping 
fan speed constant means that thrust can be rapidly mod- 
ulated without the necessity of waiting for the rotor 
speeds to change. This is especially important when 
attempting to meet the stringent thrust response requirement 
of one second between a 62% to 95% power request change. 

In addition to the steady state schedules » a number 
of limit schedules are incorporated for engine protection. 

Fuel flow is limited by maximum fan speed, core speed, and 
fan turbine inlet temperature schedules. Also fan pitch and 
nozzle area controls have maximum and minimum limit schedules. 
And finally, closed loop proportional plus integral control 
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was implemented on EPR, fan speed and inlet static pressure 
(inlet Mach number) to insure that these three outputs 
follow their respective demand schedules. 

The QCSEE control designed using the MUltivatiable 
Nyquist Array (MNA) method (described in Chapter 2) makes 
use of mnch of the control logic described above. In 
particular/ demand schedules and limit logic designed by 
the engine manufacturer are used with the MNA design as 
are the max and min limit logic. However, the MNA design 
replaces the three independently designed uncoupled servo 
loops on EPR, fan speed, and Mach number with an MNA- 
designed pre-compensator matrix^ which provides cross- 
coupling between the command schedules for optimum overall 
control performance. 

In the next section the non-linear digital simulation 
for the QCSEE engine is described. This simulation is then 
used to obtain linear models for the MNA design method. 

B. Non-Linear QCSEE Simulation and Linear Models 

During the hardware test of the UTW engine at the NASA 
Lewis Research Center, NASA developed a highly non-linear 
accurate real-time digital simvilat ion of the QCSEE engine 
at sea-leVel static conditions. This non-linear model yias 
used in extensive tests at the NASA-Ames in-flight simulator 
facility for test pilot evaluations of integrated engine air 
frame STOL configurations [17] . 


The digital simulation of the QCSEE UTW engine / shown 
in Figure 7, was based on both steady state, flow-pressure-* 
temperature relationships (maps) for each component (fan, 
compressor, combustor, etc.) and on dynamic relationships. 
These relationships are based upon continuity considerations 
for mass, momentum, and energy. 

Maps of pressure and temperature ratio as functions 
of corrected rotor speed and weight flow are used to model 
fan tip and also high pressure compressor performance. 
Compressor variable stators were assumed to be on schedule 
at all times. Pan hub performance is a polynomial function 
of fan tip performance. The inlet throat and duct are 
modeled by static pressure loss relationships. The combustor 
is represented by a heat addition and a pressure drop, with 
combustor volume dynamics neglected. Both high and low 
pressure turbines are represented by maps of flow and 
enthalpy drop which are functions of pressure ratios across 
the turbines and corrected speeds. Compressor speed is 
computed by forming a torque balance between torque generated 

by the turbine and absorbed by the high pressure compressor. 

! 

The same procedure is followed to compute fan speed. 

The complete non-linear simulation is programmed in 
Fortran and uses a modified Euler method for integration. 

It runs much faster than real-time, having a problem time to 
CPU time ratio of approximately 25 to one on a UNIVAC 1100-40 


The MNA control system design procedure requires 
linear operating point models. These models can be 
generated directly from the non-linear digital simulation 
through the use of a perturbation technique. This procedure 
(also used to obtain the PlOO linear models) [283 produces 
matrices A, B, C, and D for the standard state variable 
form model: 

X - Ax + Bu 
y = Cx + Du 


" (35) 

(36) 


where x is the system state vector of dimension eighty u 
is a third order control vector, and y is an eighth order 
output vector. 

For the QCSEB UTW engine the physical variables in- 
dicated in Table I were used to define the linear state 
models. Note that several of the variables in the y(t) 
vector are not measurable. 1 

To obtain linear models, the following procedure .was 
followed: 


1. Run the simulation to steady state 
operating point conditions. 

2 . Disconnect integrator inputs and I 

perturb each state variable Xj(t) in 

turn (j = 1, .. ., 8) by a small amount 
AXj (t) from its steady state level. 
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3r Note the cortesponding change in the 
state derivatives (AXj(t)) and outputs 

Ayj (t) . Compute elements o£ A and c 

using 

= AXj^/Ax. 


4. Compute steady state gain relationships 
between u . and x. and y^ by perturbing 

Uj in turn (j = 1, 2, 3) and allowing 

the simulation to reach steady state. 
Note the steady state relationships as 


Ax(~) = K Au(») 

. - - X. 

Ay(») = Ky Au(oo) 

5. Determine the B and D matrices as 


B = -AK. 


D = K - CK 


Following this procedure linear models were obtained 
for five power demand values, i.e., 62.5, 70, 80> 90, and 
100 percent power. The matrices obtained are indicated 
in Appendix A. 


TABLE I 


VARIABLES QCSEE^UTW LINEAR MODELS 


State Variables 

inlet duct total pressure 
Pl 3 - Pan discharge total pressure 
P^ ~ Compressor discharge pressure 
Pg - Core nozzle total pressure 
NL - Fan speed 
NH - Compressor speed 
Tg - Compressor discharge temperature 
T^ — Combustor discharge temperature 


Input' Variables 

XMV - Fuel metering valve position 

X18 - Pan nozzle area actuator position 

Theta I - Fan pitch mechanism drive motor position 


Output Variables 

^^11 ~ Engine inlet static pressure 

Pj^g - Pan discharge total pressure 

P^ - Combustor pressure 

Pg - Core nozzle total pressure 

NL - Fan rotor speed 

NH - Compressor speed 

T41C - Calculated turbine temperature 

PN - Net thrust 
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C. MNA Design 

The MNA control systems designs to be discussed in this 
section utilize the state equations of Appendix A aug- 
mented by one state variable for each integrated input variable 
This integral control action eliminates steady state offset 
and conforms to the standard QGSEE control form. Since all 
control loops utilize integral control action, the entries 
in the transfer function matrix, G(s), were modified appro- 
priately. Thus, the resulting MNA designs could be imple- 
mented directly in the non-linear simulation for design 
evaluation. 

The first MNA design for the QCSEE-UTW engine asjsumes 
nozzle area to be fixed at the full open position and fan 
pitch angle fixed at the 62.5% of full power position. With 
these two control inputs fixed at their respective values, 
the only input variable to be determined theh is fuel 
metering valve position. For this SISO condition, combustor 
exit pressure was selected as the output variable to 

be conti^olled as it is most representative of system thrust 
response. 

Prom the linear models of Appendix A, transfer functions 
augmented with integral control action for at each 

operating point are obtained. For the 62.5% condition, the 
Nyquist , inverse Nyquist, and Bode diagrams are indicated 
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in Figure 6 for a feedback gain of 0.04. The phase margin 
is 29.90 degrees with an infinite gain margin. The ciosed 
loop bandwidth at this design point is approximately 
10. radians with the feedback control 



"’.04 P- 
4 


(37) 


To examine the effect of this control at each of the 
Other four operating points (70%* 80%* 90%, and 100%) the 
composite Nichols chart of Figure 7 may be used. Here 
each curve represents the respective Nichols plot for each 
transfer function with a gain of 0.04. The 62.5% 
power curve and the 70% power curve coincide. As the power 
level increases, the corresponding curves move to the right. 

This signifies an improvement in phase margin (to 50° at 
100% power) with a corresponding decrease In loop bandwidth. 

The transient responses of the nonlinear simulation to 
the SISO design discussed above for a power change from 62% 
to 70% are indicated in Figures 8A-8D. For reference 
purposes the transient responses of the simulation using the full 
GE QCSEE control are also indicated in the figures. 


Over the power range 62.5"’80% power» the full GE control 

consists of a SISO X control and sufficient control logic 

mv ^ 


(i - e . , 
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Figure 7: Nichols Chart for SISO Design as a 
Function of Power Setting 














Figure 8C: SISO HNA, - Full C-E-0CSEE 62.5-70S 






and Theta I (fan pitch angle). When the environmental 
conditions near 80% power are realized, the GE control 
logic permits Theta I to change, at or near 90% power 
conditions is permitted to vary. It is not surprising, 
therefore, that the steady state conditions for the SISO 
MNA and full QCSEE control responses are identical over 
62.5-70% full power. 

What is surprising, however, is the improved responses 
for thrust (PN) , and fan speed (NL) of Figure 8C. Since 

total fuel consumption is the area under the Xj^^^ curve, 
it appears that the MNA control of (37) is more fuel 
efficient at this operating condition. 

To examine the impact of the constant gain control of 
(37), in light of the Nichols chart of Figure 7, a non-linear 
simulation run was made iir?r a power slam situation, i.e. , 
62.5% to 100% full power. Here the engine command setting 
of 100% full power was requested from the 62.5% operating 
point condition. The corresponding transient responses for 
the MNA-SISO and full GE controls are indicated in Figures 
9A-9D. 

Prom these responses it can be seen that fan speed (NL) 
and inlet duct pressxire significantly higher for 

the MNA-SISO design than for the full GE QCSEE control. 
However, engine thrust (FN) and total fuel consumption (area 
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Figure 9A: SISO MNA - Full 6E-QCSEE 62.5-100! 
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under X„„ curve) for the MNA SISO could be considered 

lUV 

an improvement over the corresponding full GE control 
responses. It appear Sy therefore, that both and 
Theta I (fan pitch angle are used for noise suppression 
in the full GE control and not (as stated in (15)) 

"primarily to provide rapid thrust response." Except 
for the obvious reduction in NL, the impact of Theta I 
of Figure 9C is unclear . 

For a more meaningful comparison^ the GE controls 
on nozzle area (Xj^g) and fan pitch angle (Theta I) were 
held fixed at the 62. 5% full power position. This conduit ion 
is then equivalent to the MNA SISO design and. is presented 
in Figures lOA-lOD for a transient power slam to 100% 
power. The two designs are comparable in transient and 
steady state performance. The control logic to implement 
(37) is considerably simpler than the logic in [15, 16). 

As indicated in Figures 8A-8D, good system performance 

is obtained at the 62.5% power point by a SISO control design. 

The design point at which two controls are simultaneously used 

by the GE QCSEE control occurs at the 80% power point. Using 

the linearized model corresponding to 80% power, a two input 

two output condition was configured for an MNA design 

attempt. The inputs selected were X and Theta I with 

mv 
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Figure lOA: SISO - SISC GE-0CSEE 62.5-100% 







Figure 1GB; SISO MNA - SISO 6E-0CSEE 62.5-100% 
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Figure IOC; SISO MNA “ SISO GE~0CSEE 62t5"100% 
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Here nozzle area (X^^g) was fixed at its full open position 
(62.5% power) . 

Application of the MNA design suite described in [22] 
to the corresponding transfer function matrix at 80% power 
yielded system dominance over the range 0.1 to 50 radians. 
The resulting Constant compensators are 
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Figure 14: Nichols Chart for Control Loop 2 
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of full power are indicated in Figures 15A-15D with the 
two input GE control superimposed. Nozzle area was 
held fixed at the full open position which corresponds 
to the 62,5 power setting. The peak overshoots are signi- 
ficantly lower for the MNA design with no degradation 
in performance. Also total fuel consumption (area under 
X ) is much less for the MNA design. 

For a full pov;er slam of 62.5% to 100%, the 2x2 MNA 
design is compared with the two input GE control in 
Figures 16A-16D. In every figure, the MNA design perfoirms 
very well with respect to the two input QCSEE control. 

A comparison of this 2x2 MNA design with the full GE QCSEE 
design is indicated in Figures 17A-17D. The only signi- 
ficant deviation is, as expected, in " 9^2 XMll. These 

variables are directly effected by nozzle area, 

There is, however, a significant difference in imple- 
mentation logic. This is apparent when considering the im- 
plementation of the constants in (38) -(40) versus the logic 
outlined in [15, 16] . 

To regulate inlet Mach number, a third controlled variable, 
inlet duct pressure ' is introduced. Using nozzle 

area as the third input, a three input three output design 
condition exists o Although time did not permit a full scale 
three input three output MNA design over the full operating 








Figure 15B; 2x2 MNA-Two Input QCSEE - 80-90? 
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Figure 15D: 2x2 MNA-Two Input QCSEE - 80-90% 



Figure 16A: 2x2 MNA-Two Input OCSEE - 62.5-100% 







Figure 16B; 2x2 MNA-Two Input QCSEE - 62.5-100^ 


































Figure 17C: 2x2 MNA-Full QCSEE - 62.5-100% 
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range, some analysis and simulation were performed. The 
remainder of this section is directed to this condition. 

Since the GE nozzle area control Was designed using 
SISO methods, the SISO MNA and the 2x2 MNA designs were 
independently augmented with the GE nozzle area control. 

The first case represented in Figures 18A-18D consists 
of the SISO MNA with the GE control While Theta I 
is fixed at. the 62.5 power position. Note the rapid 
thrust response (same as in Figure llC) and the respective 
decreases in P ^2 XMll from their counterparts in Figures 
9A and 9D, Execept for fan speed, the Theta I curve for the 
full GE control has little influence on the remaining 
variables. 

Auginenting the GE nozzle area control to the 2x2 MNA 
control design, results in the transient performance 
indicated in Figures 19A-19D. As expected, the dynamics 
are unchanged from the previous case except for the improved 
performance of P ^2 ^^NA control performs very 

well with the GE nozzle area control for noise suppression. 

I 

Preliminary design using MNA concepts for the three 
input three output case was performed using the 90% power 
point. The input vector u(t) and output vector y (t) were 



















Figure 18D: SIS0 MNA with GF-MB Control - Full OCSEE - 62.5-100% 
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Figure 19A: 2x2 MNA with GE Nozzle Area - Full 6E Control 62.5-100% 
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Figure 19D; 2x2 MNA with GE Nozzle Area - Full 6F Control 62.5-100% 



chosen as 






u(t) = 

^18 

y(t) e 

^12 


Theta I 




where is nozzle area and Inlet duct pressure. 

For this case, system dominance was easily obtained 
for each linearized operating point in Appendix A. The 
CAD procedure of [18 , 22] resulted in the following 
compensators when 0.1 < w < 40 radians at the 90% power 
point 


K 


. 2032E-4 

.0300 

.0289 


-.01856 
“7.9359 
103 « 67 


0,00157 

“.02023 

“.001195 


71.455, 0.4523} 

P = DiagdO.O, 0.1, 20.0} 

The corresponding DNA, image, and Bode diagrams for this 
operating point are presented in Figures 20-22, 
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The MNA designs discussed above demonstrate the utility 
of the design concept in achieving good turhofan engine 
performance. The full potential of the Bode and Nichols 
concepts for multivariable plants were not fully utilized 
in this design application, however » as dynamic compensators 
were not required. In each case, appropriately determined 
constant compensators provided adequate closed loop system 
performance and feedback gain scheduling was not needed, 

This desirable condition is due principally to the well be- 
haved and easily controlled dynamic system structure. 
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APPENDIX A 

LINEAR MODELS FOR 
THE QCSEE ENGINE 

1. 62.5% power level - Model includes actuator states. 

2. . 70% power without actuator states. j 

3. 80% power without actuator states. 

4. 90% power without actuator states. 

5. 100% power without actuator states. 
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